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Abstract. The classical theory of strain in material continua is reviewed and 
generalized to space-time. Strain is attributed to "external" (matter /energy fields) 
and intrinsic sources fixing the global symmetry of the universe (defects in the 
continuum). A Lagrangian for space-time is worked out, adding to the usual Hilbert 
term an "elastic" contribution from rntrrnsic strain. This approach is equivalent to a 
peculiar tensor field, which is indeed part of the metric tensor The theory gives a 
configuration of space-time accounting both for the initial inflation and for the late 
acceleration. Considering also the contribution from matter the theory is used to fit 
the luminosity data of type la supemovae, giving satisfactory results. Finally the 
Newtonian limit of the theory is obtained. 
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1. Introduction 

There are in cosmology a limited but important number of facts which do not fit in 
the classical general relativistic (GR) frame of the standard theory. Among others we 
recall the current accelerated expansion of the universe, as discovered at the end of 
1998 observing the redshift dependence of the luminosity of type la supernovae |lll|2|], 
the extreme uniformity of the temperature of the Cosmic Microwave Background 
(CMB) [3} the anomalous rotation curves of spiral galaxies lHHZl, the anomalous 
mass-luminosity relation in galaxy clusters llSj and the galaxy power spectra [9]. If, 
on one side, the number of facts to be accounted for is limited, on the other, the 
number of theories put forth to explain them is extremely high. All these theories 
may however be grouped under a limited number of headings. One is GR plus 
appropriate fields. We may include in this category the inflaton, accounting for 
the homogeneity of the CMB (lOl HH |12l, as well as quintom QSl [HI, k-essence 
|fT5l[T6B and vector inflation |[T7|. Then we find dark energy in a number of variants: 
cosmological constant IMIIlUlQllZll, quintessence |l22l HI/ phantom energy |l24ll25|. 
holographic dark energy | |26ll27ll28| , triads |29| and k-essence again [30J. Together 
with the idea of dark energy we may quote an, in principle, different ingredient of the 
universe which is dark matterfST]; now the competition opens among many possible 
types of non-baryonic matter and specific candidates: Hot Dark Matter [32J, which 
actually seems not to be the right stuff since it does not allow for structure formation; 
Cold Dark Matter ||33l : Warm Dark Matter p4|; neutrinos, gravitinos, photinos, 
axions, sterile neutrinos. Weakly Interacting Massive Particles (WIMPS), Lightest 
Supersymmetric Particles (LSP). By the way the most successful cosmological theory 
of the moment is the so called Lambda-Cold Dark Matter theory (ACDM) ||35l 

Another compartment is the one of extended, alternative, modified theories of 
gravity (all changes are with respect to standard GR). In this group we find both 
metric and Palatini versions of / {R) theories ||36l |37| and Modified Newtonian 
Dynamics (MOND) [SSHSSHQI and its covariant formulation Ell. 



All these possibilities remain within the realm of classical theories, but of course 
we have also a host of attempts to introduce quantum aspects in a more or less 
fundamental way. So we find various string cosmologies ||42|, brane world ||43| , Loop 
Quantum Cosmology | l44ll45| , and similar approaches. 

Everything can be good or not according to the capability of one or another 
idea to account for a consistent enough number of facts, which, as said, are rather 
few. The impression one has, however, is that most theories rely much more on the 
mathematics of the model rather than on physics. In some cases we have more or less 
explicit tricks introduced in order to produce the right equations. The wording is in 
any case physical, but the entities that are called in often have rather exotic properties, 
in the sense that they are completely unusual, i.e. not corresponding to anything else 
we know. So, considering the various entities that are introduced in the universe 
in order to account for the accelerated expansion we have "fluids" which do not 
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dilute even though the available room is increasing (cosmological constant), violating 
the strong energy condition (quintessence), such that when compressed they expand 
(phantom energy). Mathematically everything is OK but someone may feel uneasy 
with all that belonging to the actual world. 

In order to find a route in the wild forest of conjectures one needs a criterion 
and one could be to stick, as far as possible, to the physical properties one already 
knows and has learned to describe. Of course there is no reason why the universe on 
the largest scale should possess properties that look like the ones we are used to at 
everyday's scale, but we have some hints telling us that this approach could be not 
completely absurd. On one side we know that the world becomes richer at the tiniest 
scales than at our usual scale: wave functions and the whole machinery of quantum 
mechanics are not needed to describe the world at the human scale. On another side 
we know that some properties of complex systems are scale invariant. On this more 
or less inspiring basis we have tried to draw on a known methodology and a known 
theory in order to describe space-time as such. In GR space-time is indeed described 
as a real continuum which interacts with something else, being this "something else" 
what we call matter/ energy; the interaction is expressed by the non-linear Einstein 
equations whereby the non-linearity implies that space-time, in a sense, interacts also 
with itself. Following on we easily see that many techniques used for space-time in 
GR are very similar to the ones applied to describe material continua in the theory of 
elasticity: in the end, in both cases, everything is dominated by geometry. How far 
this analogy can go? How the basic tools of the elasticity theory need be modified in 
order to fit the behaviour of a four-dimensional manifold with Lorentzian signature? 
The idea is not at all new and in a sense has distant roots even in the old investigations 
on the nature of the "luminiferous ether" that haunted physics until the beginning of 
the XXth century. Of course the "ether" was thought of as something filling space 
and moving or vibrating in an absolute time, the mathematical tools and concepts 
were the ones of the XIX century and the aim of the research was the explanation of 
the propagation of light. However, as we know, after getting rid of the name "ether" 
Einstein remarked (1920) that "...space is endowed with physical qualities; in this sense, 
therefore, there exists an ether..." |j46l, and in today's quantum field theories the ether 
is named "the vacuum" and has indeed important physical properties. So what we 
are about to do is to explore the idea of space-time as a material continuum using 
the formalism and general framework of GR and the leading idea of the material 
continua we know. The central concept will be the strain of the four-dimensional 
manifold which is related to the metric tensor and is in turn produced by some cause. 
The origin of the strain can be either "external" or internal and in our case, according 
to the universally adopted dualistic approach, "external" will mean matter/ energy. 
The possibility of internal or spontaneous strain states is well known in the theory 
of elasticity and plasticity and is connected with the presence of structural defects in 
the "texture" of the manifold. So the other concept we shall use in our theory will be 
the one of space-time defect or Cosmic Defect (CD) which we use also to name the 
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theory. 



We have already applied this approach heuristically to the interpretation of the 
late acceleration of the universe, obtaining encouraging results | |47ll48| . Here we are 
giving a new formulation of the theory, casting it in a systematic and general way 

2. Review and reformulation of the metric properties of elastic continua 

For our purposes we shall consider an N— dimensional material continuum and 
assume that in its reference state it is perfectly homogeneous and isotropic so that 
we may identify it with an N— dimensional Euclidean manifold (reference manifold) 



If we now imagine to apply some action (in three dimensions we would speak 
of forces) either localized or diffused, we expect the continuum to be brought to a 
deformed (strained) state, where in general it, as a manifold, will acquire curvature: 
this will be our "natural" manifold |50| . 

Let us then imagine we have two copies of the system, one in the reference 
state, the other in the natural one, and that they are both embedded in an N + 
1 - dimensional simply-connected flat Mtnkowskian manifold [J the additional 
dimension will be time-like. We consider a situation in which it is possible to put 
in one to one correspondence points on the two manifolds. This operation may be 
formally thought of as a type II gauge transformation [52ll53| . 

Considering the embedding N + 1 - manifold we may introduce a "displacement 
vector field" u which connects pairs of points in the two embedded N -dimensional 
manifolds. If r is the N + 1- vector localizing a point in the reference manifold, the 
u vector runs from that position to a different one r' in the natural manifold. The 
situation is shown on fig.([l]). X's are the (Cartesian) coordinates of the embedding 
manifold; the coordinates on the reference manifold are ^'s; the ones of the natural 
manifold are x's. We assume that the functional dependences leading from one 
coordinate system to the other are all sufficiently smooth for the rest of our reasoning 
to be valid. Of course the different dimensionality implies the existence of at least one 
constraint for each of the two sub-manifolds allowing for the dimensional reduction. 
In general we assume = (X^,..., X]v+«) and (on the reference sub-manifold) 
X; = hi (Xi, X]v) for z = N + 1, N + n ; the same holds for the natural manifold: 
= Woe 0^1, ■■•/ ^N+w) and X/ = qi (Xi, Xjv) for z = N + 1, N + n. If all functions 
are differentiable and invertible one can also write: 



:j: Actually for general natural manifolds it can happen that more than 1 extra-dimensions are needed 
in order to build a flat embedding manifold [ ,51 J- In particular, it is known that an embedding of a non- 
vacuum solution in five dimensions is minimal while an embedding of a non-flat vacuum solution in 
six dimensions is minimal. 



ins. 




(1) 



A tensor theory of sf ace-time as a strained material continuum 



5 




Figure 1. Embedding of the reference and the natural manifolds in an N + 1 - 
dimensional flat manifold. Xq are the Cartesian coordinates in the embedding 
manifold; are the coordinates on the reference (sub)-manifold (N - dimensional), 
Xy are the coordinates in the natural (generally curved) (sub)-manifold. u represents 
the displacement vector from points of the reference manifold to points of the natural 
manifold. 



On these bases we may write, with reference to the embedding manifold: 

/ (X) = r (X) + u (X) (2) 

The functions / and w together with the constraints allow to cast ^ in terms 
of the reference {^} or natural {x} coordinates. A choice one can always do is 
to numerically identify ^'s and x's for corresponding places in the two manifolds, 
not forgetting that this choice can be practically useful but has per se no physical 
meaning. 

Everything is non-trivial only in the case of u being a non-trivial field: rigid 
translations are uninteresting and amount to a simple coordinate change in the same 
manifold. 

The next step we may think to do is to compare the distances between pairs of 
corresponding points in the unstretched and stretched situations. This comparison is 
meaningful only if made in one and a single manifold, i.e. the embedding one. Using 
unprimed symbols for the reference manifold and primed symbols for the natural 
one we may write 



V2_^ A^ia^^ih a,b = l,...,N + n 



df = r/aidX^dX^ 
dl'^ = t]^i,dX'^dX" 

In both cases the metric tensor is the one pertaining to the flat embedding 
manifold. Introducing the constraints defining the two sub-manifolds of interest and 
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using the coordinates adopted for each of them, the equations above become 

Saturation of the Latin indices leads, on the first row, to the N -dimensional flat 
metric tensor of the reference manifold (either Euclidean or Minkowskian); however 
it is important to stress that the additional factors 9^"/9x'' do not correspond to a 
coordinate change on the same manifold, so that the final j^^ tensor is not in general 
a metric tensor on the natural manifold. Viceversa on the second line, the constraint 
to stay on the natural manifold leads in general to a curved N - dimensional metric 
(with Lorentzian signature). Using ^ and ([T]) we see that it is: 

g}iv (x) = + 2e^v (x) (3) 
where ^^y is the strain tensor of the natural manifold, given by 

_ 1 / du" du^ du"du^\ 

Looking at dU) we easily see that the strain we have defined transforms as a 
genuine tensor on the natural four-dimensional manifold. 
If one could find a coordinate transformation such that © could be written as 

then the strained and the unstrained situations would be diffeomorphic to each other 
(j^v would coincide with g^y) and it would be impossible to perceive the deformation 
from within the manifold: no intrinsic curvature. In fact the integrability condition 
for ^ is De Saint Venant's: 

R^s = 0' (6) 
where R'^r^g are the components of the Riemann tensor. This is indeed the case of 
elastic deformations, which have an exogenous origin: they are due to the application 
of external forces and the strain is brought to zero whenever those forces are removed 
(absence of plastic deformations). In practice a smooth and continuous u field leads 
to an integrable (O. If we are interested in intrinsic deformations (the ones that 
can be sensed from within) we must study singular u (and e) fields, the singularity 
being represented by some discontinuity in u and/ or its derivatives. Now, a singular 
displacement field in a continuum means that the medium contains what is formally 
defined as one defect (or more), according to the definition given by Volterra in 1907, 
while studying elastic and plastic deformations [54] . 

The singularity in u reflects of course also in the strain tensor and in general we 
shall write |55l l56l |57| the elementary deformation (from ((Z) in the reference frame) 
as a non integrable one form: 

dr'^ = coUr\ 
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so that the deformed line element will be 

ds^ = giiv (x) dx^dx^ = rji^^colcOpd^^d^P 

Now, our target is space-time and we know that interesting situations there 
imply condition © to be violated so that we are led to the conclusion that any 
non-trivial space-time should contain at least one defect in the sense recalled above. 
This is a different way to the singularity theorems by Hawking and Penrose [58 J. A 
warning that is appropriate to issue at this point is that the defect we are considering 
here should not be confused with the topological ones often appearing in cosmology. 
There the topological defects are the residual of the phase transition which gave 
origin to matter during inflation. Of course "our" defects also have topological 
properties, however their nature is completely different, as we have seen. 



3. Lagrangians for space-time 

After the geometric considerations developed in the previous section we are left 
with the problem of finding the functional dependence of u (or to say better, e) on 
the coordinates. In practice this means that we need to introduce an appropriate 
Lagrangian for our space-time manifold. 

Let us start from the typical Einstein-Hilbert Lagrangian for a defectless 
manifold 

Leh = Ry^. (7) 

which is build from the simplest scalar obtainable from the curvature tensor. We of 
course expect that the Lagrangian we are looking for reduces to © when the defect 
disappears. An apparently reasonable approach is to add to ^ an "elastic" potential 
term: this would be consistent with the description given so far. 
A typical expression for the elastic potential energy is: 

We = \cr^yt^' (8) 

where the usual meaning of cr-^iv is to be the components of the stress tensor of the 
material. The stress is naturally dependent on the strain and viceversa. In the 
standard elasticity theory this dependence is mostly dealt with assuming a linear 
relation (Hooke's law) [49J. In the case of space-time we have a priori no reason to 
say that it is so also, but we shall assume linearity and see what happens. We write: 

where Caii^v are the components of the elastic modulus tensor and are supposed to 
be independent from the e's. As a consequence © becomes: 

W, = ^C,,,^£-/^£?- (9) 

We have of course to do with a scalar density and we may study the situation 
in a locally flat manifold (tangent space). Assuming that the medium, at least in the 
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unperturbed condition, is perfectly homogeneous and isotropic the form assumed by 
the elastic modulus tensor depends on two parameters only 11491 and is 

Coi^^S = ^Vct^VjS + V {Voc^V^S + VasV^^) (10) 

The parameters A and }i are known as the Lame coefficients; their value is a property 
of the continuum under consideration. In the standard elasticity theory instead of 
the Tj's one has Kronecker deltas (at least for Cartesian coordinates), i.e. the Euclides 
metric tensor; here we introduce the Minkowski tensor, because we are dealing with 
space-time. In the case of space-time the global isotropy condition in four dimensions 
has to be taken with caution, because of the light cones, however we shall assume it 
holds. Considering ((TO)) and ^ we arrive at 



2 

being e = the trace of the strain tensor. In the globally curved manifold the 
corresponding Lagrangian density is: 

he = QAe2 + ^e.^e-f^^ ^\ (11) 

Of course indices are raised and lowered by means of the global metric tensor. 

The complete action integral for the natural manifold (no external forces: in 
practice no matter, in our case) will be: 

j [r + \xz^ + ^e.^e'^^^ sj\^\dx^ (12) 

The structure of ((12)) recalls the classical form where a kinetic and a potential 
terms appear. Here the role of kinetic term is played by which contains derivatives 
of the strain tensor. 



3.1. "Elastic" Einstein equations 

The treatment we have given of the behaviour of the "elastic-style" space-time is not 
different in the form from the introduction of some "material" contribution to the 
Lagrangian. Now this contribution is expressed by (|TT|) with 
1 

Sfiv = 2 {S}iv - ifiv) ■ (13) 

One has to be careful in dealing with ]piy. As already said, this is a tensor on 
the natural manifold, but is no metric at all for that manifold. The j^^ tensor can in 
practice be identified with the metric of the local tangent four-dimensional frame, 
comoving with the cosmic flow of the given universe. At different cosmic times the 
various local j^^ tensors are related to each other via boosts based on the expansion 
rate. 
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In principle from ((TTj) and using ((131) , varying the action integral with respect to 
g^y (or equivalently e^y which is proportional to the non trivial part of the metric 
tensor), we may obtain generalized Einstein equations in the form: 

G^v = Tepiv + kT^v (14) 

where all that is neither in the Einstein tensor nor in the matter energy-momentum 
tensor can be interpreted as an effective "elastic" energy momentum tensor Tg^v. Since 
Teiiy, as well as T^y, is obtained varying a true scalar (the integrand of ((T2|) ) with 
respect to a true tensor, it is also a good tensor, retaining all properties of tensors. In 
particular no coordinate choice can bring Tg^y to zero, unless it is identically zero. 

In vacuo the final matter term is absent, however now it is in general R ^ 0, 
provided some internal cause, like a defect, is there. The tensor Te^y, though being 
partially built from the metric itself, plays the role of an additional source together 
with the proper matter term. In vacuo, for instance, we see that the Bianchi identities 
applied to the Einstein tensor, imply that the "elastic" energy-momentum tensor is 
conserved. When matter is also present, the conservation condition applies to the 
sum Tepiy + KTpiy and, in general, the possibility of a transfer of energy between the 
matter and the strain term is given; this is not different from the mechanisms which 
can subtract energy from material systems pumping it into gravitational waves in 
GR, but for the fact that now the energy of the wave has a new gauge independent 
interpretation. To the whole source (T„v and Te^y together) we may apply the energy 
conditions often considered in GR; if the matter source is thought to be an ordinary 
one the conditions are separately satisfied for it, so that we obtain constraints for the 
Lame coefficients of space-time. 

Our procedure starts from an assumption of homogeneity and isotropy for 
space-time. The latter constraint does not come from the specific type of universe one 
wishes to implement, but from the properties of the reference unstrained manifold. 
In fact, in an unstrained manifold, where neither defects nor material sources are 
present, no event and no direction is better of any other, so the system is unavoidably 
homogeneous and isotropic. When considering a peculiar global configuration, 
induced by a defect or an arbitrary matter distribution, the corresponding locally 
anisotropic strain of the real manifold can of course induce in turn some anisotropy 
(and inhomogeneity) in the elastic parameters, however, as it is the case for ordinary 
three-dimensional continua, the induced anisotropy in the properties of the "stuff" 
may usually be considered as second order with respect to the strain. So, as far as the 
theory is assumed to be linear, the homogeneity and isotropy of the elastic modulus 
tensor can be held. 

4. Specific symmetries: a Robertson Walker space-time 

The purpose of the previous section was to find out an appropriate Lagrangian for 
space-time and it explicitly depends on the strain tensor which, in turn, depends on 
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the symmetry. In the present work we are interested in describing the universe as a 
whole so we now focus on Robertson Walker (RW) symmetry and derive the strain 
tensor in such a symmetry by means of eq. (|T3l) . Willing to use the displacement 
vector field (cfr. eq.H]) we need an embedding flat manifold. From a geometrical 
point of view the RW symmetry may be interpreted as a point symmetry so that 
a natural choice is to have polar coordinates around the center of symmetry. Even 
though the system is four-dimensional, it is instructive, as a first step, to describe 
it in two dimensions: one radial distance, which for us will be cosmic time, and an 
"angle". So, in the spirit of the previous sections, our reference manifold is a plane 
(with Lorentzian signature) and the natural manifold is a curved bi-dimensional 
surface (again with Lorentzian signature). Both surfaces are embedded in a flat 
three-dimensional manifold and the natural choice of coordinates for it will be the 
cylindrical ones. With no loss of generality we may assume that the reference 
axis passes through the center of symmetry of the natural manifold. The global 
configuration is reproduced in figlZl The T coordinate is the radial coordinate in 




r 



Figure 2. A curved surface with a central symmetry is embedded in a three- 
dimensional manifold. The reference frame is a plane and the global coordinates are 
cylindrical. 

the natural manifold. The reference manifold is localized by means of the constraint 
z = constant; the constraint for the natural frame is z = / (r) being / a non-linear 
function, otherwise one would have a cone, i.e. a flat sub-manifold. The global 
coordinates are r, a, z; in the reference frame they are r, a coinciding with the first 
pair of global coordinates; in the natural frame the coordinates are t and a, with 
a coinciding with the corresponding global coordinate. Using the flatness of the 
embedding manifold we see that 

dT^ = dr^ + dz^ = (l + f^) (15) 
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Of course /' is a shorthand notation for df / dr. If / is a regular function (with the 
possible exception of the origin) from (|T5|) we can work out r [t). 

It is easy to write down the distance between two nearby points on the reference 
frame directly adopting the appropriate signature: 

ds^ = dr^ - r^dcr^ + dz^ = dr^ - r^dcr^ (16) 

The distance between the corresponding points in the natural frame is 

= dT^ - r^da^ = + f'^) - ^^^^^ (1^) 

We can read out the metric tensor on the curved manifold from (|T7l). both in the 
r, a and the t, a coordinates, and, just to recover the formal correspondence with the 
RW notation, we may identify r with a{r). 

The natural coordinates for our curved manifold are T and a, so, in order to 
proceed, we also need to cast (|T6l) in terms of these coordinates; then: 

dT^ 
1+/ 

The difference between ((T7|) and ((181) is 



ds^ = ^r-^ - (18) 



- ds^ = f'^dr^ = ' dT^ (19) 



l+/'2 



Interpreting (|T9|) on the light of sect. (O we read out the only non-zero element of 
the strain tensor of the natural frame and using the coordinates thereupon (use (flSll ) 



/'2 _ 1 



fl2 



2(l+/'2) 2 

With a little change in the notation we put a = dr / dr; the reason is to conform to the 
standard notation for a RW universe. 



4.1. 4 + 1-dimensional embedding 

The above treatment, when applied to the full four-dimensional manifold we use 
to describe the universe, actually corresponds to a negatively curved space (space 
curvature constant k = —1), as can be seen by noting that the choice z = constant 
gives a flat Lorentzian manifold only if dcr^ = dx^ + sinh^;^ (^d9^ + sin^ 9d(p^), which 
corresponds to a three-dimensional negatively curved space. Of course we would 
like to analyze the null space curvature case since cosmological observations point in 
that direction. Let us then look a little more in detail to the geometry of our manifolds 
and to the meaning of the cr variable we used in the previous section. Actually 
various embedding strategies have been adopted for similar purposes 1*5911511. We 
shall stay with our "cylindrical" symmetry approach and write a five-dimensional 
flat line element in the form: 

dsl = dr^ - dal + dz^ (20) 
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Here dcfQ is the three-dimensional space line element in the null space curvature case: 

dal = dr^ + {dO'^ + sin^ ed(p^^ ; 

z is the fifth coordinate of the embedding space-time. With this coordinate choice, 
if we put z = constant we recover a 4-dimensional flat sub-manifold, while RW 
symmetry with A: = is obtained with the following transformation, as shown in 
Ell EH: 



a[T) r 
T(r,T) 
z(r, t). 

Here a{f) is the scale factor in the four-dimensional RW metric, and the explicit 
transformation for t and z are 

dt = ^fl + ^ + flr^^ rfT/2 + fl(T)rrfr 

dz= {^ — \— flr^^ dx 17. — a{T)rdr. 

The final purpose is to work out the strain tensor that, using (|T3l) , turns out to only 
have spatial components: 

{0 if pi = V = 1 

if ^ = v = 2,3,4 
if V. 

We can now compute the trace and the second order scalar that appear in the elastic 
Lagrangian ((TT|) : 

£a/3g ^ - e - 2 — ^2 — 3 

In the case of non-null spatial curvature the dependence of the strain tensor on the 
scale factor is different according to the sign of the curvature parameter. In practice 
the embedding strategy requires different coordinates and different transformations 
for each k value. Recalling sec. we see that the strain tensor both for positive and 
negative spatial curvatures has only one component different from zero, namely the 
time-time component: 

erT = \{l + ka^), k = ±l. 

In the case of a negative space curvature this has been explicitly worked out in 
our 2+1-dimensionaI example. For the k = 1 case, however, the reference four- 
dimensional flat manifold has to be Euclidean, even though the natural manifold 
has a Lorentzian signature. This is due to the fact that the 3-dimensional spatial sub- 
manifold is a 2-sphere. 
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4.2. The space-time behaviour 



Let us now maintain a RW symmetry and study the case of an empty space-time in 
which space is also flat {k = 0) as apparently it is for the actual universe. The action 



integral (|T2l) becomes 

S = 



■ 2 

aa + a 



a + Le 



with 



dr, 



8^- 



(21) 



2>2 



where B = A + |^ is the bulk modulus of the continuum. Here we are using a point 
Lagrangian and we shall derive the equations of motion varying the action in eq.(f2T|). 
We can remark that, in the absence of a defect, a RW space-time would in general not 
be a solution of the Lagrange equations obtained from ((T2)l: the symmetry (and the 
defect) are a priori conditions thus forcing the action integral ((2T]) . 

Since second order derivatives in the Lagrangian appear linearly we can get rid 
of them by means of an integration by parts so that the effective Lagrangian density 
becomes 

L = 6aa^ + Lg (22) 

We can work out the energy function, W = — L, which is, by construction, a 
conserved quantity: 



W = Gad 



■2 



2>2 



Solving (|23|) for a we have 



6a 16 



2^2 



(23) 



(24) 



We see that for a ^ one has a ^ oo. For a ^ co a diverges also. 
The expansion rate in (|24)) has a minimum, which means that at the beginning 
the expansion is decelerated, then it becomes accelerated. 

The trend of ((24l) appears in fig.© for arbitrary positive values of both parameters. 
In order to recover General Relativity in the absence of defects, the energy function 
must be set to zero so that from now on we assume W = 0. 



5. Expansion of the universe 

The main fact we want to account for is the accelerated expansion of the universe, so 
in the present section we will deal with the data evidencing this phenomenon. We of 
course shall proceed according to our approach, but in the following we would also 
like to show how the same formulae can be read and interpreted in more traditional 
ways. 
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Figure 3. Expansion rate of a Robertson Walker space-time versus the scale factor, 
according to the theory. The space is assumed to be flat. Arbitrary positive values of 
the parameters have been used. 



5.1. Dark fluid interpretation 



Let us put together the first integral in (|24l) and the second order evolution equation. 



deduced from the point-Lagrangian in (|22|): 



„2 ^ 3 



16 



-Tda — a 



a'' 

16 a 



2^(3fl2 + l), 



(25) 
(26) 



H is the Hubble parameter. Interpreting the r.h.s. of the equations as representing a 
fluid component, we may read out the corresponding density and pressure: 



16 fl4 
3 '-^ 



a 

,2 . 

16" ~fl4 



a^-1) 



(27) 



(28) 



The state parameter, i.e. w = p/p, would clearly depend on time. Since we are 
interested here in late cosmology, let us derive the condition for the acceleration to 
occur. An accelerating phase, a > 0, requires |0 + 3p < 0. In our model it turns out to 
be 

so acceleration sets in when a^ > 1, or, in terms of the redshift z = ao/a — 1, when 
z < ao — 1. The parameter is the present scale factor and its value depends on the 
model and the observation. 

In particular, if we write down the equation for the "elastic" state parameter 



P 

w = — 



1 3fl2 + 1 
3 
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we can easily see that the behaviour of the elastic potential tracks radiation, curvature 
and cosmological constant as a increases, passing from w = 1 /3 in the a ^ limit 
to IV = —1 for ^ oo. In a 3+1 view we may think that, close to the cosmic defect, 
a release of "elastic" energy in form of radiation (primordial gravitational waves) 
dominates; afterwards this radiation is progressively converted into the equivalent 
of a "dark energy" driving an accelerated expansion, just as a cosmological constant 
would do. 



5.2. Type la supernovae luminosity 

The most direct evidence for the acceleration of the expansion comes from the 
luminosity data from the type la supernovae. In order to test the theory on the Snia 
data we must include in our analysis the presence of matter. This will be done as 
usual introducing in the Lagrangian a matter term minimally coupled to geometry. 
It is 



The coupling constant k is 16nG/c^ and Wi is determined by the equation of 
state of the ith component; pio is the mass/energy density of the ith component in 
the comoving frame, the index refers to present day values. The luminosity is 
commonly expressed in terms of the distance modulus and the redshift parameter 
appearing in the scale factor through a = aq/ (1 + z). It is ||60| 



2\2 3{1+Wi 



ar 



^_M = 25 + 51og((l + z) /^'^) (31) 

For (|3T1) to hold, distances have to be measured in Mpc. 

The simplest is to restrict to dust and radiation, for which iv is respectively and 
1/3; re-organizing the constants, the Hubble parameter is 



'l + z)2\^ 8k 



3 1- ^ o M +—{1+Z)^p,n0+Pr0{l + Z)] (32) 



^\ V~^^J 3b' 
Using (|32|) , eq. (|3T|) gives the explicit form of the distance modulus. 



5.3. Fitting the data 

Once (|32l) has been introduced into (HD we may use an optimization procedure in 
order to fit the luminosity data from type la supernovae. The implied range of values 
of z (<^ 2) is sufficiently small to assume that the radiation contribution is negligible 

(Pro ^ 0). 
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The distance modulus (|3T]) can be written in terms of three optimization 
parameters and the new integration variable ^ = (1 + z) /aq / as in the following 



m - M = m + 5 log f (1 + z) 



3 (^/2_ 1)2 + ^^/3 



(33) 



where 



m = 25+ |logl6- |logB + 51ogflo 



3B "O 



(34) 



We have compared the results in (|33|) with the data of 307 Snia's from the 
Supernova Cosmology Project union survey [611 • The best fit is shown on figlH 




Figure 4. Fit of the luminosity data from 307 type la supemovae obtained applying 
the CD theory. Three optimization parameters have been used. The reduced of the 
fit is 1.017. 



The optimal parameters values are: 
m = 46 ± 1 

flo = 1.97 ±0.04 (35) 
xp = 10.1 ±0.3 

Considering the different sensitivity with respect to changes in m, which is outside 
the logarithm in (|33l) , and in arid ip , which are in the logarithm, the optimization 
has been performed in two steps. In the first step the optimization routine ^621 has 
been run with all three parameters giving the actual value for m and a first estimate 
of flo and xp (actually 1.98 and 10.18) with a big uncertainty (up to 100% in the case of 
ip). In the second step m has been fixed to its already found value and the routine has 
been used again with the two remaining free parameters, thus yielding a far better 
uncertainty. The final reduced of the fit is 1.017. 
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We compare our result with what can be obtained using the ACDM theory. In 
this case the distance modulus, expressed in terms of two free parameters, is B48l 

dz' 

^ + Slog (1 + z) + Slog / ^^^^=^^^^= (36) 

•'^ /q^ (1 + z') V 1 - n„, 



m - M 



Using the same two-steps optimization process the final is 1.019, so we see that 
the CD fit is better. 

Let us now deduce some cosmological parameters from our best fit values by 



using (|34|) and (|35|) . Eliminating B between the first and the second equation of (|34)) 
we find a relation among all parameters we have used for the fit, and po- 

m = 25 + ^log-^. (37) 

The actual value of the mass/ energy density in the universe is not easy to assess, but 
is usually thought to be somewhere between 10~^^ and 10~^^ kg/ m'^. 
By using eq. (|37|) we can derive the estimated mass density from our fit, and it is 
Pq ^ 3.4 X 10~^^kg/m^. This is but a rough estimate because the Snia data do not 
allow for more accurate results. Actually the uncertainty domain for pQ is in the order 
of 100%. Another estimate we can do concerns the value of the Hubble constant Hq. 
From ((32)), neglecting the contribution from radiation, we have: 



, , 1 (l-alf 1 

Ho = Cv^y^^-^ + g (38) 

Introducing the numerical values we found for the parameter^ we have 

H„ = 64±35-^. 

s X Mpc 

Once more the big uncertainty comes mainly from the uncertainties in the luminosity 
data of Snia's. 

Until now we have avoided the explicit use of the B parameter, however this 
parameter, in our theory, has the simple physical meaning of bulk modulus of space- 
time, so let us compute it. From the first equation of the (|34|) we obtain: 

B = (3 ± 2) X 10"^Mpc"2 = (3 ± 2) X 10~^^m-^. (39) 

Looking back at equation ((24)) we find for B the constraint B > 0, which is indeed 
satisfied by (l39b . 



6. Newtonian limit of the theory 

After having found a good correspondence between theory and data at the level of 
the SnIa luminosity dependence on redshift we should also verify that a correct weak 
energy (Newtonian) limit exists. To that purpose we start with some general remarks. 

§ Remember that for the use in the magnitude calculations all distances must be expressed in Mpc. 
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Since in practice our theory simply additively introduces a peculiar source term into 
the Lagrangian of space time, whenever this new term is sent to zero we of course 
recover plain GR with all its features and local limits. It would not be so only if the 
additional term (|TT|) were somehow singular, which is not the case. In fact, excluding 
the very cosmic defect, he can continuously go to zero at any place, together with the 
local strain. 

We think this could be enough, however let us explicitly verify what the weak 
field limit is. In order to make this check, let us consider a spherically symmetric, 
stationary physical system. We know that the general line element for this problem 
is 

ds^ = f (r) dr^ - h (r) dr^ - r^dO^ - r^ sin^ Odcp^ (40) 

where Schwarzschild coordinates are used. 

It is easy to read our strain tensor (its non-zero elements) out of the metric in (|40)) 
comparing it with a Minkowski line element in spherical coordinates. We get: 

T :s («) 

Equipped with (|4T]) we are able to explicitly work out (|TT|) for our problem, then 
we write the corresponding Lagrangian for empty space-time (no matter, except at 
the source of the local spherical symmetry). Modulo an unessential sinO factor and 
after eliminating a second order derivative of / by a by part integration, the explicit 
formula is (primes mean derivatives with respect to r): 

^ _ Ifh'r+lfh^-lfh A (/-l)(fe-l) ^.2 



From (|42|) we obtain the Euler-Lagrange equations for / and h: 



fh-rf'-f 



rh' + h^-h 



16fh 



A 



^2h + lf-5)f-Afh^ + {f-hy 



+2fi [ (ih^ + 2h-3^f + {l- If) /z^] I 



16/2 



A 



(4/ -3)h^ + {2h - iff - 2fh 



(43) 



(44) 



The left hand side of eq.s (|43ti44|) corresponds to the equations for the 
Schwarzschild problem, whilst the right hand side can be considered to be "small" 
as far as it is 



r^^,r^\ << 1 



(45) 



< 
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Under the latter assumption we may look for solutions like: 

/ = /o + /i (46) 
h = ho + hi 

where 

is the Schwarzschild solution (m = GM/c^) and 

fiM-K^ (47) 
Under the ansatz (|46l) and condition (|47| eq.s (|43ti44[) may be solved to first order 
in A, }i. The result is: 

- ;2^Hfy + (22^ - 6^2) In (r - 2m) - f mr (l - 2^) ) 
hi = — -C ( t^Kn, , -m^r- 4m3 In (r - 2m) ) 

+ , ^ f v^T^ - ^^^^ - - lO'^^ In (r - 2m) 

Here A and C are integration constants. A is used to remove constant contributions 
from fi. 

Up to this moment we have not used the assumption that the gravitational field 
is weak. Let us introduce this condition now, linearizing the results in m /r. The 
result is 

/i ~ ^ - l^mr 
hi ~ 7 + 4^ - \}imr 

As for the integration constant C, it contributes a redefinition of the mass of the 
spherical source and a short range contribution to hi; fixing C = we finally have: 

-l-^i-l^mr 

= 1 + 2^- \}imr ^ ' 

The solutions (|49|) are acceptable, as already said, as far as condition (|45l) holds. 
Now condition (|45)) can be updated to 

1 

r << — 
m^ 

If we consider the example of a Sun-like star, where it is m ~ 10'^ m, and use the 
result found for B in the previous section (reasonably it is also fi ^ B) we obtain 

r « lO^^m IQ^Mpc 10^^ lighty ears 

As we see all deviations of our theory from the standard GR, at the Solar System 
or galaxy clusters level, are absolutely negligible. The only effect is at the cosmic 
level. 
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7. Conclusion 

We have exploited the existing analogy between the theory of elasticity and GR and 
this approach has given good fruits, however, recalling the open questions we posed 
in the Introduction, we remark here that an analogy is not an identity; we are not 
allowed to mechanically transpose one theory on top of the other. In GR one properly 
looks for "static" solutions in four dimensions: various possible configurations of the 
full four-dimensional universe are studied, but there is no evolution because there 
is no evolution parameter out of the manifold. Time is part of the manifold and the 
dynamic term in the Lagrangian is the scalar curvature which contains second order 
derivatives of the Lagrangian coordinates, i.e. of the elements of the metric tensor, 
with respect to an arbitrary set of Gaussian coordinates. In the case of elasticity 
the three-dimensional manifold has Euclidean signature and time is the absolute 
Newtonian time; the dynamics is expressed by time derivatives. The difference is 
paramount. 

The role of the defect(s) in our theory deserves also some additional comments. 
There are indeed, in the three-dimensional theory of continua, situations where one 
has to do with continuous distributions of (micro)defects, rather than with localized 
defects. This happens mainly as a consequence of a plastic deformation: an initial 
stress in the material is eased (for instance by reheating) and produces a distributed 
defects field. The final result of this process is indeed the disappearance of the 
internal strain and a permanent plastic deformation of the manifold. This situation 
does not correspond to what we want to describe, since for us the strain represents 
the gravitational field, or the non-trivial part of the metric tensor, and we want it 
to stay there and disappear only when its causes are removed (no plastic shear). 
This is the reason why our attention is focused on a localized defect responsible 
for a global spontaneous strain state and for the related symmetry. An appealing 
possibility could be to have a plurality of localized defects which would give rise 
to some more complicated strain pattern; one could be tempted to identify these 
defects with matter. We have not dared, for the moment, to pursue this idea, so in 
our theory, at present, matter is treated as usually it is, appearing in the Lagrangian as 
an additional independent term coupled with the (strained) geometry via the metric 
tensor. 

In the above conceptual background we have combined the classical GR 
approach with the description of space-time by means of the linear theory of 
elasticity, preserving general covariance and all the features of GR. The expansion 
of the universe is consequently described in terms of a strained four-dimensional 
continuum (space-time), whose strain has partly an intrinsic origin due to the 
presence of a cosmic defect, partly depends on extrinsic sources, i.e. matter fields. 
The cosmic defect fixes the global symmetry of the universe, then the general features 
of what, in our (3 + l)-split view, is the cosmic expansion. This approach can 
be applied to any kind of universe with any global symmetry depending on the 
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possible defects, just as it happens in three dimensional solid materials. We have then 
assumed the basic "stuff", i.e. space-time, to be locally homogeneous and isotropic, as 
a consequence of the global homogeneity and isotropy of the flat unstrained reference 
manifold and of the linearity of the theory. Working out the global configuration of 
space-time for a RW symmetry (the typical symmetry assumed to hold, on a cosmic 
scale, for our universe) we have found that it naturally includes an initial extremely 
rapid expansion with a steeply decreasing expansion rate, followed by acceleration. 
We have also verified that the inclusion of matter in the form of fluid(s) preserving the 
global symmetry does not modify the general structure of the expansion. The theory 
depends on three parameters, which are the present scale factor of the universe ao, 
the bulk modulus of space-time B, and the present day matter/ energy density of the 
universe pQ. Using these three quantities as optimization parameters we have fitted 
the luminosity data of type la supemovae. The result has been good, and the value 
obtained for jOq is consistent with the current estimates for barionic matter without a 
need for more matter, but the uncertainty due to the accuracy of the luminosity data 
is very high. According to the theory no further dark energy is needed; however, if 
we wish, we may read our elastic contribution as a dark energy fluid, whose density 
and pressure have been explicitly written, it would however be rather difficult to find 
a reasonable physical interpretation for the properties of this peculiar fluid. 

A final remark concerns the signature of our manifolds. It always is Lorentzian 
in the natural manifold which we want to represent the actual universe. As for the 
reference manifold it can either be Euclidean or Minkowskian and the embedding 
strategy can easily produce one or the other of them. It is sufficient to assume the 
embedding flat higher dimensional manifold to be Minkowskian. Whenever then 
the reference submanifold is a space-like hyperplane (time-like normal vector) its 
geometry is naturally Euclidean; viceversa, choosing as a reference submanifold a 
time-like hyperplane, it will turn out to be Minkowskian. 

Although the theory has been applied to the cosmic scale, we have also verified 
that it has a correct Newtonian limit and, with the value of the parameters obtained 
from the Snia's fit, it is indistinguishable from GR at the Solar system as well as at 
the galaxy clusters scale. 
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